Strategies for enhancing quantum entanglement by local photon subtraction 
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Subtracting photons from a two-mode squeezed state is a well-known method to increase entan- 
glement. We analyse different strategies of local photon subtraction from a two-mode squeezed 
state in terms of entanglement gain and success probability. We develop a general framework that 
incorporates imperfections and losses in all stages of the process: before, during, and after subtrac- 
tion. By combining all three effects into a single efficiency parameter, we provide analytical and 
numerical results for subtraction strategies using photon- number-resolving and threshold detectors. 
We compare the entanglement gain afforded by symmetric and asymmetric subtraction scenarios 
across the two modes. For a given amount of loss, we identify an optimised set of parameters, such 
as initial squeezing and subtraction beam splitter transmissivity, that maximise the entanglement 
gain rate. We identify regimes for which asymmetric subtraction of different Fock states on the two 
modes outperforms symmetric strategies. In the lossless limit, subtracting a single photon from one 
mode always produces the highest entanglement gain rate. In the lossy case, the optimal strategy 
depends strongly on the losses on each mode individually, such that there is no general optimal 
strategy. Rather, taking losses on each mode as the only input parameters, we can identify the 
optimal subtraction strategy and required beam splitter transmissivities and initial squeezing pa- 
rameter. Finally, we discuss the implications of our results for the distillation of continuous- variable 
quantum entanglement. 

PACS numbers: 03.67.Bg, 42.50.Ex, 03.67.Hk, 03.67.Pp 



I. INTRODUCTION 



Efficient distribution of entanglement between distant 
parties is fundamental to most quantum communication 
protocols. However, entanglement is fragile and suf- 
fers from decoherence, which is detrimental to the per- 
formance of any communication protocol upon which 
it relies. The ability to increase the entanglement be- 
tween communicating parties is therefore vital, and fur- 
ther, practical considerations dictate that this should be 
achieved through only local operations and classical com- 
munication (LOCC). One solution — entanglement dis- 
tillation — extracts a small ensemble of more strongly 
entangled states from a larger ensemble of weakly en- 
tangled states [ ]. In the discrete variable regime, en- 
tanglement distillation has been achieved using photonic 
qubits [ ]. In the continuous- variable (CV) regime, the 
situation is more involved. Most common CV states and 
operations are Gaussian in nature. However, there exists 
a no-go theorem which states that one cannot distill en- 
tanglement from Gaussian states by Gaussian operations 
alone [ ]. Gaussian operations are those with Hamilto- 
nians which are (at most) quadratic in the ladder opera- 
tors a, eft, comprising the basic tools of quantum optics 
including beam splitters, phase shifters, squeezers and 
homodyne detection. A number of protocols to increase 
entanglement in CV systems have been proposed [4-7], 
elements of which have been implemented [8-12]. 
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Photon subtraction was initially proposed by Opatrny 
et al. to increase the efficacy of a teleportation proto- 
col [4]. Since then, several studies have looked at photon 
subtraction in more detail. Cochrane et al. investigated 
subtracting and detecting n photons simultaneously from 
the modes of a two-mode squeezed state [13]. Olivares 
et al. studied the use of on-off (single-photon thresh- 
old) detectors to measure at least one photon subtracted 
from both modes coincidentally, again in terms of the 
improvement of a teleportation protocol [14]. Kitagawa 
et al. provided a detailed numerical analysis of two-mode 
subtraction by on-off detectors in terms of the explicit 
change in entanglement and compared this with the op- 
erational measures used previously in the literature [15]. 
This work was built on by Zhang and van Loock [16] 
in which analytical results for perfect symmetric sub- 
traction using photon number resolving detectors and 
on-off detectors were derived. More recently, Navarrete- 
Benlloch et al. obtained similar results and quantified 
the non-Gaussianity of the operations [17]. 

Until now, efforts have been focussed on symmetric 
subtraction, that is, identical measurement results fol- 
lowing subtraction of photons from both modes of an 
entangled state. In this paper, we ask for the best en- 
tanglement enhancement strategy based on local photon 
subtraction from each mode independently. We inves- 
tigate if there is a difference in entanglement gain or 
success probability when detecting different numbers of 
photons subtracted from each mode, which we refer to as 
asymmetric subtraction. Perfect symmetry between the 
modes is impossible to achieve in practice. In this paper 
we account for asymmetric losses and detector efficien- 
cies, which enables our results to be applied directly to 
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realistic experiments. 

The article is organised as follows: in Sec. II, we define 
the figures of merit by which protocols that increase en- 
tanglement may be measured. In Sec. Ill, we derive the 
state evolution of photon subtraction from a two-mode 
squeezed state (TMSS) and we provide an analytical form 
for the probabilities corresponding to different detection 
strategies. We use the evolved state to analyse the gain 
in entanglement by different subtraction strategies: using 
photon-number resolving detectors (PNRDs) in Sec. IV, 
and threshold detectors (such as avalanche photodiodes 
— APDs) in Sec. V. For both types of detector, we nu- 
merically analyse symmetric and asymmetric subtraction 
in the presence of loss occurring before, during, and after 
subtraction, and show how these different parameters af- 
fect the success of the protocol in Sec. IV D. To aid the 
reader, the locations of the principal results and the main 
figures in which they are depicted are summarised in the 
table below. We also list the main conclusions drawn 
from our work in Sec. VI. 
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II. ENTANGLEMENT, GAIN AND RATE 

Entanglement in this system can be captured conve- 
niently by the positive partial transpose (PPT) criterion. 
If the partial transpose p TA of a density matrix p has neg- 
ative eigenvalues, then p must be entangled [18, 19]. The 
sum of the absolute values of the negative eigenvalues of 
p TA is defined as the negativity N (p). The entanglement 
measure we use in this paper is the log-negativity [20], 
defined 



E N (p) =log 2 [l + 2JV(p)] =log 2 



(1) 



where ||X| 



Tr 



denotes the trace norm of X. 

The maximally entangled CV state, for a fixed en- 
ergy, is the two- mode squeezed state (TMSS), as gener- 
ated during parametric down-conversion (PDC) or equiv- 
alently by interfering two single-mode squeezed vacua in 
phase at a 50:50 beam splitter [21]. This state can be 
written in the Fock basis as 



I^tmss) = Vl - \ 2 ^\ n \n,n) A A> 



(2) 



n=0 



^l-X>Y^—a) n a'^\Q,Q) AA , , (3) 
z — ' n! 



n=0 



which describes n pairs of photons in modes A and A' 
for a given squeezing parameter A £ [0, 1) and we define 

\n,m) A A' = \n)A ® \m) A f - 



The density matrix of this state is 

oo oo 

Ptmss = (1 - A 2 ) Y, E A"A TO |n) A (m|®|n)^(m| , (4) 

n=0 m=0 

the partial transpose (with respect to mode A) of which 
is 

oo oo 

P™ss = (1 " A 2 ) E E * n nrn)A(n\®\n) A (m\ . (5) 

n=0 m=0 

Taking the trace norm yields 
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(1 - A 2 ) 
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n=0 



Tr 



X> m |m)<r 



.m=0 



\2 ' 



(6) 



(1-A) 

therefore the TMSS defined in Eqn. (2) has log-negativity 

1 + A 



En (ptmss) = log 2 



f - A 



(7) 



which is the benchmark against which changes in entan- 
glement will be measured. 

From this initial state Ptmss, a subtraction step s pro- 
duces a state p 5 with entanglement En (ps)- We de- 
fine the gain in entanglement G (p 5 ) as the difference 
between entanglement after distillation and that of the 
initial state, normalised to the initial entanglement, i.e. 



G(p 5 



En (p s 



En (ptmss) 



1 , 



(8) 



such that G (p s ) > if the subtraction step increases 
entanglement. 

The probability of an entanglement enhancing step 
must be less than unity, since entanglement cannot be 
increased on average under LOCC [2 ]. It is calculated 

by 



P(p,)=Tr\s(p)] , 



(9) 



where 5 (p) is the unnormalised density matrix following 
a subtraction operation. The final density matrix p 5 is 
found by normalising $ (p) thus 



Hp) 

Tr [*(/>)] 



(10) 



The aim of an entanglement enhancement protocol is 
to increase entanglement from some initial value. It is 
therefore desirable not only for subtraction to produce a 
high gain, but also to do so at a high rate, i.e. that the 
probability of a subtraction event is as high as possible. 
In general these conditions cannot be met independently, 
therefore there exists an optimum to be found based on 
the parameters of a given implementation. We therefore 
define the entanglement gain rate (frequently shortened 
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FIG. 1: A two-mode squeezed state (TMSS) in a photon- 
subtraction setup. A TMSS initially occupies modes A, A' . 
Following the array of beam splitters (with transmissivity 
Ti,T/), we seek the entanglement in modes Si,S[ following 
detection in modes 63, 63. Losses are modelled by photons 
in modes before (5*2,6*2), during (54,64) and after (65,65) 
detection. 



to "rate" in the remainder of the paper) T (p s ) as the 
gain afforded by a distillation step G (p s ) multiplied by 
its likelihood 



T( P ,) = P(p,)G(p,) 



(11) 



By post-selecting on successful subtraction events, the 
entanglement of this sub-ensemble can be increased. The 
increase will depend on the parameters of the subtraction 
step employed, and can be maximised to determine the 
values of those parameters which yield the highest rate. 



III. THE PROTOCOL 

To effect subtraction and account for loss, we consider 
an array of 8 beam splitters, each with transmissivity 
Ti and reflectivity R { = ^l-T 2 , acting on the TMSS 
as shown in Fig. 1. The state initially occupies modes 
A, A! . The beam splitters T 2 , T' 2 effect subtraction, while 
the other beam splitters Ti^T^^T^ model losses before 
subtraction, during detection and after subtraction, re- 
spectively (and similarly for the primed counterparts). 
As such, the state vector describing the combined state of 
the input modes may be written |^i n ) = I^tmss) ® |0)® 8 . 
The evolution of the input modes is governed by the uni- 
tary transformation U U' . The elements of [/, U' de- 
pend on the arrangement of beam splitters and how they 
couple their respective modes. For the array of beam 
splitters shown in Fig. 1, U may be written as 
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(12) 

and similarly for U' with all symbols replaced by their 
primed counterparts. Since the only non- vacuum input 
modes are A, A', the output modes depend only on the 



first column of each component unitary U : U' . Indeed, 
we may write the output state vector in the Fock basis 
as 



n=0 ' \m=l 



E<m<C) |0> 



(13) 



where ai^ m is the (1, m)-th element of the unitary U and 
1S the creation operator for the mode m. We are in- 
terested in the entanglement between the modes S± , S[ , 
occupied by e, e' photons respectively, dependent on de- 
tecting <i, d! photons in modes S3, S3, respectively. The 
modes S2 , 64 , S5 , S' 2 , S4 , S5 contain Z, Z' photons lost to 
the environment. These loss modes can be combined, 
simplifying our problem into the following circuit dia- 
gram 
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The photons in two input modes of the TMSS are divided 
into three output modes each: entangled, detected and 
lost, and we assign parameters a, a', /?, j3' and 7, 7' to be 
the fraction of photons in each mode, respectively. These 
may be written in terms of the components of the unitary 
U thus 



'1,1 
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tItItI , 

2 7~>2rr-i2 



F = < 3 = TtRp£ , 
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(15a) 
(15b) 
(15c) 



m=2,4,5 



where the coefficients satisfy a 2 + ft 2 -f 7 2 = 1, as defined 
by the unitarity condition (and similarly for their primed 
counterparts). Substituting Eqns. (15) into Eqn. (13), 
yields 

|*out> = y/l - A 2 J2 A "E c e,dA>,d',i'\ e > e '>K d')\l, V) , 



n=0 



(16) 



where the summation is over S = {e, e', <i, d', Z, V s.t. e 
d + I = e' + d' + V = n}, and where we have defined 




C e ,d,l = 



with the multinomial coefficient 



n 
e, d, Z 



(17) 
(and sim- 



ilarly for the primed quantities). On tracing out the loss 
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and detected modes, we obtain the (mixed) entangled 
state p across modes Si : S[ as 



CO ^max,^ max Z max 

P=( 1 ~ x2 ) £ S H c «,»,<i,<i',;,;'|e)s 1 (e|®|e , )s;(e , | 

n,n=0 

(18) 

where d max = min (n, h) — / max = min (n,n), and 
Cn,h,d,d' are subject to the constraints e + d + / = 
e' + d' +V = n, e + d + / = e' + d' + /' = n, such that 



C"n,n,d,d' ,1,1' 



\ n+n n\h\a e+i a' e ' +s ' /3 2d f3' 2d ' ' 1 2l 1 ' 21 ' 
d\d'\l\l'\Vd^W\ 



(19) 

As the mode structure is now clear, we will suppress the 
mode labels Si, S[ in further discussions for compactness. 



A. Subtraction 

Subtraction is effected by changing the limits of the 
summation over d, d! in Eqn. (18), leading to an unnor- 
malised subtracted state s(p). The exact form of this 
state depends on the subtraction strategy employed and 
is derived in detail in the relevant sections below. To 



calculate the entanglement we seek the sum of the nega- 
tive eigenvalues of the partial transpose of the normalised 
state. Using Eqns. (10) and (1), we therefore seek 



E n (p s ) = log 2 



*(P) 



Tr[ S (p)] 



(20) 



where 5 (p) Tsi is the partially transposed unnormalised 
state with respect to mode Si. Following the analysis 
in Refs. [15, 16], we can write s (p) Tsi in block diagonal 
form, and we change variables such that 

n z + Z+d, n j+i+d, I' i+j-K+l + d-d' , (21) 

where the indices z, j denote the rows and columns of the 
K th block matrix of dimension K + 1. Thus s (p) Tsi is 
written explicitly as 

oo K 

S (p) Ts i = (l-A 2 ) C\f\m®\K-i){K-j\ , 

K=0i,j=0 

(22) 

with coefficients 



J ~ d ht t> h Hl + i + j + d-K- d')\d\d'\ v /i\j\(K-j)\(K-i)\ 



where Iq = max {0, K + d! — d 



the matrices 



C 



i,3 



— j}. By inspection 
are sym- 

,K;j=0,...,K 



amounts to removing the sums over d, d! in Eqn. (23). 



metric, and they are also persymmetric (and therefore 
centrosymmetric) when the primed parameters are equal 
to their unprimed counterparts and the range of the sum- 
mations over d and d! are equal [ ]. Eqn. (23) is central 
to the derivation of all subsequent analytic results. 

The coefficients lead to a normalised state if 

d, d! are summed over entirely, i.e. the limits on the 
summation over d, d! are t = t f = 0, t max = t' max = 
oo respectively. The (unnormalised) state 5 (p) fol- 
lowing a subtraction event is found by placing limits 
on the d, d! summation, yielding coefficients Cl[^ = 



C^*P (£, £max 5 Cax)) dependent on the strategy em- 
ployed. For threshold detectors that click on receipt 
of a minimum number of photons, the sum runs from 
the threshold value t to t max = oo (for instance, single- 
photon avalanche detectors have t = 1). For photon- 
number-resolving-detectors (PNRDs) operating within 
their resolution regime, the summation disappears since 

^max = t- 

IV. SUBTRACTION USING PNRDS 



A. Probability of photon subtraction 



The probability P (t, t') = P (p s ) of detecting t = t max , 
t' = t' max photons, respectively, is given by the trace of 
the unnormalised state 5 (p) after a subtraction event. 
Detecting t, t' subtracted photons by PNRDs projects 
onto modes Si , S[ the unnormalised state 



oo K 

B(p)= (l-A 2 ) C\f\i) E {j\®\K-i)E>{K-j\ , 

K=0i,j=0 

(24) 

where the coefficients = (£,£',£,£') are iden- 

tical to those in Eqn. (23) for t max = t, t' max = t' and 
therefore the summation over d is dropped. This has 
probability 



We will first briefly consider the case when sub- 
tracted photons are measured using perfect PNRDs. This 
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FIG. 2: Lossless symmetric subtraction of a single photon using perfect photon- number resolving detectors. Entanglement 
gain (a), distillation probability (b) and rate (c) are shown for A G [0, 1) and a 2 G [0, 1). 



P(M')=Tr[ 5 (p)] 



(l - A 2 ) \ 2t f3' 2t f3 2t (a 2 + 7 2 ) 
[l-(a 2 +7 2 )(a /2 + 7 /2 )A 2 ] 



t'-t 

P (*'-*,o) 

t'+i r t 



l+(a 2 +7 2 ) (a /2 +7 /2 )A 2 
1- (a 2 +7 2 )(a /2 + 7 /2 ) A 2 



(25) 



C { *\ This yields 



where P^' 5 ^ [z] denotes the n th -order Jacobi polynomial. while the denominator traces over the diagonal elements 
We have assumed, without loss of generality, that t <t'\ 
owing to the symmetry of the problem a corresponding 
equation for t > t' can be found by swapping all of the 
primed values with their unprimed counterparts. This 
expression is valid if t, t' are known, i.e. in the case where 
detection resolves photon number. This may be extended 
to threshold detectors by summing over t, t' from the 
threshold t,t' (typically unity) to WjCax- We deal 
with this scenario in greater detail in Sec. V. 



J- Tr [j(*>C<*> 



B. Lossless symmetric subtraction 

Perfect symmetric photon subtraction detects an equal 
number of photons in both modes with 100% efficiency. 
This is equivalent to setting 7 = 7' = 0, whereby I = V = 
and a = T2 , the transmissivity of the subtraction beam 
splitter. Symmetry implies t = t', a = a' and /3 = 
implying K = i + j, which simplifies Eqn. (23) to 



C ( K ) = x K+2t a 2K /3 4t (i + ty.(j + t y. 

1,3 WyJi\j\{K-i)\{K-3)\ 



°K,i+j 



(26) 



where the Kronecker delta function $K,i+j means that all 
off anti-diagonal elements are zero. In this case, each of 
the (K + 1) x (K + 1) blocks C (K) are both symmet- 
ric, = Cj*\ and persymmetric c\ K ^ = C ( £\ K _ j , 
which can be exploited [16, 24] to compute the log- 
negativity using 



E N (p) = log 2 



(l-A 2 )E^ Tr[J^)C^)] 
P(t,t) 



(27) 



where 3^ K > is the exchange matrix and we include ex- 
plicitly the normalisation factor in the denominator. The 
numerator traces over the antidiagonal elements C\ I ^_ i 



K=0 i=0 

it 



m\(K-i)\ 



A/3 2 \ 



1 



^2 ' 



,1-Aa 2 / (l-\a 2 Y 
and from Eqn. (25), the probability simplifies to 



1 - A 2 A 2t /3 4t 

P(t,t) = ± — i-rrPt 

K ' (1 - A 2 a 4 ) t+1 



l+(Xa 2 Y 
1 - (Aa 2 ) 2 



(28) 



(29) 



where P n [z] is the n th Legendre polynomial, a special 

case of the Jacobi polynomial P^ 1 ^ [z] found in Eqn. (25) 

given by Pn°'°\z] = P n [z] . This is an analytic, closed- 
form expression of the results in [15, 16], and yields an 
entanglement of 



E N (p s ) = log 2 




1 - (Aa 2 ) 2 



(30) 

Fig. 2 depicts the result for varying squeezing parame- 
ter A and subtraction coefficient a 2 . The behaviour of 
the entanglement gain rate in Fig. 2(c) shows that there 
exists an optimum value that provides the best entangle- 
ment yield per trial. High gain is less likely, such that 
the rate is peaked at particular values of A opt = 0.66 and 
a opt = 0-83. On average, these parameters produce the 
highest gain in entanglement per trial. 

In general, the parameter a 2 is freely tunable when 
performing the protocol, whereas A is restricted by the 
maximum squeezing available. Therefore it is useful to 
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obtain an expression for the optimal subtraction parame- 
ter a 2 pt in terms of the initial squeezing A that maximises 
the entanglement gain rate T (p s ). These values are ob- 
tained from Fig. 2(c), and are plotted as a function of A 
in Fig. 3(a). For t = 1, the fit to these data is a second- 
order polynomial of the form 



a 2 opt = 0.238 (A - l) 2 + 0.576 (A - 1) + 1 . (31) 



Thus, given that in the lossless case a = T2, we have the 
recipe to dial up the most effective subtraction rate to 
maximise the gain in entanglement per trial. 

It has been shown previously that all things being 
equal, subtracting more photons increases the gain after 
subtraction [ ]. However, by optimizing a 2 with respect 
to the rate, this gain is only marginally higher with in- 
creasing <i, as shown in Fig. 3(c). Indeed, at A > 0.6, 
there is no advantage to subtracting more photons. As 
expected, even with the optimal a 2 , the probability of 
entanglement gain decreases significantly with increasing 
t, shown in Fig. 3(b). These effects combine to yield a 
lower entanglement gain rate as the number of subtracted 
photons is increased, as shown in Fig. 3(d). 



C. Lossless asymmetric subtraction 



We now consider the case of perfect asymmetric sub- 
traction, i.e. detecting different numbers of photons t ^ t' 
in each subtracted mode. When counting resources in 
terms of total photons subtracted, asymmetric subtrac- 
tion allows us to double our space to include odd photon 
numbers. This additional degree of freedom lifts the de- 
generacy in the coefficients, i.e. the primed quantities 
may take different values from their unprimed counter- 
parts. This doubles the number of parameters over which 
the protocol can be optimised. 



Starting again from Eqn. (23), we detect £, t' photons in 
each mode with unit efficiency. Since we are neglecting 
the effects of loss in the system, 7 = j' = and as 
such the only non-zero contribution to the summation 
over I is the I = term. Photon number resolution is 
maintained by setting t max = t, t' max — t' in Eqn. (23) 
as before, however they need not be equal as in the case 
above. In this scheme the relation 7 = 7' = ensures 
that all elements of Q^ K \ except those that satisfy i-\-j = 
K + t r — t, are zero. We may thus express the elements 
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of C^> as 

C (K) =(1 _ x 2 )x i +j+ 2t a i +ja ,2K-i- jp 2tp,2t> ^2) 



tWly/i\j\(K-i)\(K-j)l 



H+j,K+t'-t • 



The matrices only have elements along one of their 
skew diagonals, shifted from the main skew diagonal by 
t f — t. This allows us to define an antidiagonal subma- 

trix B( K ), where b( k ) is a (K + l) x (K + l) matrix 
and where K = K — t' + t, which contains the ele- 
ments on its main skew diagonal. Without loss of gen- 
erality, we are able to choose t < t f ', whereby the ele- 
ments of B(^) are B^ = C$,_ tJ+t ,_ t . This is both 
symmetric and persymmetric, so the above approach 
can again be used in computing the negativity. When 

the elements of interest are 



evaluating Tr 



b(*)j(*) 



B 



i+t'-t,K- 



,(*) 



where 



2t' a K+2{t'-t) n ,K R 2t R ,2t' 



(i + t')\(K -i + t')\ 



(33) 



tW\y/{% + t'- t)\(K -i + t'- t)\(K - i)W. 
Eqns. (33,25) together lead to the relations 

, „ (1 - A 2 ) A 2t '/3 /2t '^ 2t a 2 (*'~*) 

2 n l2\2\ t '+ 1 



(1 - a 2 a ,2 A 2 ) 



(t'-t,o) 



1 + (aa'X) 
1 - (aa'X) 2 



(34) 



and 



]T IV fe( R ) j( R )j =(1 - A 2 )(A/3') 2t '/? 2t a 2 ( t '- t ) 



if=0 



(i + f)l 



\2 = 



y/i\(i + t' - t)\t\t'\ 

(35/ 



Unfortunately, no analytical expression could be found 
for the sum of series in brackets. However, by considering 
the ratios between successive terms in the series, it is 
clear that it converges over the relevant range < \aa' < 
1. This gives an exact, but not closed, form for the log- 
negativity of the state 



En(Ps) =log 2 



(l-aV 2 A 2 ) £/+1 



I 1 — a 2 a f 



2 A 2 
2 A 2 



21og 2 ^(Aa^ 



.2 = 



y/i\(i + t' -t)\t\t'\ 



(36) 



1. Comparing symmetric and asymmetric subtraction 

As with the symmetric case earlier, we can directly 
compare the entanglement gain, probability and rate for 
asymmetric subtraction. For illustration purposes, we 
compare different ways of subtracting up to two photons 
from each mode. The results are shown in Fig. 4. 

It is clear from Fig. 4(a) that the symmetric cases, 
(t,tf) = (1,1) and (2,2), produce more gain than their 
asymmetric counterparts for a fixed number of subtracted 
photons. From Fig. 4(b), there exists a regime of A > 0.7 
for which the probability of subtracting (2,0) photons is 
greater than the (1,1) case, therefore the gain rate, shown 
in Fig. 4(c) is correspondingly higher in this regime. This 
is because the (2,0) subtraction event is more likely than 
the (1,1) for high values of A. However, it is clear that 
in general the asymmetric (1,0) case produces the most 
entanglement gain per trial, once again due to the high 
likelihood of subtracting 1 photon as opposed to 2 pho- 
tons in any combination. 



D. Losses and imperfect detection 

In any realistic scenario, where the process of entan- 
glement distillation and entanglement enhancement will 
be most essential, the detectors are imperfect and evolv- 
ing quantum states suffer losses. Losses are accounted 
for by setting 7, 7' > 0, and the summation over I (which 
denotes the number of photons lost) includes contribu- 
tions from I > 1 which are non-zero. Losses are modelled 
by beam splitters before, during and after subtraction, 
as shown in Fig. 1. The effect of loss on the gain re- 
gion (for the symmetric cases where one (blue), two (red) 
and three (green) photons are detected in each mode) is 
shown in Fig. 5. As the combined loss 7 2 increase, the 
range of initial squeezing A for which entanglement can 
be enhanced by local photon subtraction decreases. Fur- 
thermore, detecting higher numbers of photons is more 
loss-tolerant; since a smaller a 2 is required to subtract 
more photons, the restriction on 7 2 < 0.5 is reduced. In- 
deed, we can define 7m ax as ^ ne maximum loss for which 
entanglement still increases on subtraction. The depen- 
dence of the maximum losses 7 ma x = 1/t + 1 comes from 
the gain as a function of loss, setting to zero and solving 
for t the difference of Eqns. (39) and (7). 

The constituents of 7, 7' are determined by Eqn. (15c), 
where we define efficiencies before, during and after de- 
tection as T 2 , T 2 , and T 2 respectively. From the ex- 
pressions of 7,7', it is clear that each of these efficiencies 
contribute differently to losses. The scaling of each form 
of loss is shown in Fig. 6, for one (a) and two (b) pho- 
tons subtracted symmetrically. While the scaling with 
each form of loss does not change, it is clear that the 
protocol becomes more loss-tolerant when subtracting 
more photons, consistent with Fig. 5. It is apparent from 
Fig. 6 that the effect of detector efficiency, parameterised 
by T3, is almost negligible. This can be interpreted by 
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AAA 
(a) (b) (c) 

FIG. 4: Lossless symmetric and asymmetric subtraction with PNRDs. Using a, a optimised over entanglement gain rate for 
different combinations of subtracted photons (£,£' < 2), we calculate (a) the achievable gain, (b) the success probability and (c) 
entanglement gain rate as a function of A. 




A 0.75^0.25 



FIG. 5: Lossy symmetric subtraction with PNRDs. Region 
in a, 7, A parameter space for which gain is positive, for one 
(blue), two (red) and three (green) photons are subtracted. 



considering the subtraction detector as something of a 
post-selector: it may not often click, but when it does 
one can be fairly certain that entanglement has been in- 
creased. Losses after subtraction may be mitigated in a 
similar way, by using a loss-tolerant entanglement detec- 
tion strategy, such as post-selection or a loss-tolerant en- 
tanglement witness [ ]. Loss before subtraction, which 
may be modelled as a mixing of the state due to non unit 
channel transmission, is unavoidable and its contribution 
is significant. 



1. Entanglement gain and rate under loss 

To study the enhancement of entanglement in the pres- 
ence of losses and imperfection, we again start from 
Eqn. (23), which in the symmetric detection PNRD case 
for t = t' , is given by 



< K ) - yi+j+2{t+l) n 2K R U^2{i+j-K+2l) 
l=lo 

(i + Z + £)!(j + Z + £)! 



(37) 



+ i + K)\t\ 2 ^i\j\(K-j)\(K-i)\ ' 



where Zo = max (0, K — i — j). Since this matrix is 
symmetric and per symmetric, we can calculate the en- 
tanglement following [16]. Defining x = 1 — a 2 A and 
y = [a 2 + 7 2 ) A allows us to write 



Tr 



J(K) C (K)- 



(1 - A 2 ) (3 4t X 2t 



> 2 - 7 4 A 2 ) t+1 
and the success probability is given by 



x l + 7 4 A 2 



7 4 A 2 J ' 



1 - A 2 p At \ 2t 

P(t,t) = ~ -/--LI Pt 

V ; [1 - y 2 \ 2 ] t+1 
This yields a log-negativity of 



1 



E N {p) = log 2 



i-y 2 



7 4 A 2 



1 



x 2 - 7 4 A 2 



(38) 



i+y 2 

l-j/2 



(39) 



These results reduce to the lossless symmetric detection 
case given by Eqns. (28, 29) when 7 = 7' = 0. 



2. Optimising parameters under loss 

Subtraction is a general strategy for increasing entan- 
glement. However, it is interesting to investigate which 
two-mode squeezed states, parameterised by A, are most 
improved by this scheme. Furthermore, using Eqns. (38) 
and (39), we can incorporate loss as an additional input 
parameter. For a given implementation, the experimen- 
talist has free choice of the reflectivity of the subtraction 
beam splitter (proportional to a). Therefore a logical 
question to ask would be: given losses of a certain level, 
what is the optimum range of A for which subtraction 
yields the highest entanglement gain rate, and to what 
value should the subtraction beam splitter reflectivity be 
set in order to achieve this? 

Figure 7(a) shows the pairs of A and a 2 values re- 
quired to produce the highest entanglement gain rate in 
the presence of losses. Finding an empirical fit to these 
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A 



(a) (b) 

FIG. 6: Lossy symmetric subtraction with PNRDs. Region in loss-parameter space for which gain is positive for a range of A 
and subtraction beam splitter transmissivity T 2 = 0.83. The contribution to the overall loss parameter 7 is from the individual 
constituent efficiencies before, during and after subtraction, , T|,Tf , respectively. The scaling with each type of loss is the 
same for (a) 1 and (b) 2 photons subtracted, however, in accordance with Fig. (5), the range is larger when more photons are 
subtracted. 




(c) (d) 

FIG. 7: Lossy symmetric subtraction with PNRDs. (a): squeezing A and subtraction a 2 parameters, optimised to produce the 
highest entanglement gain rate at a given loss; and the resulting (c) gain, (b) probability and (d) rate resulting from these 
parameters. 

data one obtains the following formulae for setting the optimal A and a 2 : 

«o P t (7 2 ) = e- 381 (- 2 + 01 ) - O.67 2 + 0.8 , (40) 
A op t (7 2 ) = e- 107 - 1 ^ - 1 ) 2 + e -2-8(-r 2 +o.i) _ Q 2 _ (41) 
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The gain and rate resulting from these parameters are 
plotted as functions of the loss parameter 7 in Figs. 7(c) 
and 7(b) respectively. One can immediately see that 
there is an upper threshold on loss, above which entan- 
glement cannot be increased. 



be a multiplexed array of n APDs with a detection event 
triggered off n-fold coincidences. 



A. Probability of photon subtraction using APDs 



3. Symmetric subtraction, asymmetric loss 

By setting t' = t in Eqn. (25), one may calculate the 
probability of symmetric subtraction allowing for differ- 
ent parameters a, /?, 7 on each mode. However, lifting the 
degeneracy of the parameters across the different modes 
means the submatrices C^ K ^ are no longer centrosym- 
metric, and therefore the entanglement is not analytically 
tractable. Instead we calculate numerically the entangle- 
ment from Eqn. (23). Fig. 8 shows the gain, probability 
and rate for non-degenerate losses when subtracting a 
single photon from each mode. The results in Fig. 7 are 
simply line-outs for 7 = 7'. 



E. Losses in asymmetric subtraction 

When subtracting photons asymmetrically, the effect 
of loss on each mode may not be equivalent. As the sim- 
plest example of asymmetric subtraction, we consider the 
(t,t f ) = (1,0) case, whereby one photon is subtracted 
from one mode, and the other mode is left unchanged. 
As shown in Fig. 4(c), this case produces the highest en- 
tanglement gain rate for all A. We calculate numerically 
the optimal parameters a^ pt , A opt that give the largest 
rate of gain, as shown in Fig. 9. Note that, since we do 
not subtract from the second mode, the optimal a' 2 is, 
trivially, 1 — 7 /2 for all 7 2 . 

The resulting gain, probability and rate are shown in 
Figs. 10(a), (b) and (c) respectively. As is to be expected, 
the effect of loss is not symmetric in this case. The gain 
region is bounded more sharply by losses in the unsub- 
tracted mode, whereas the losses in the subtracted mode 
follow the scaling of the symmetric case, Fig. 7(a). 



The four measurement outcomes to consider when us- 
ing APDs are as follows: both detectors clicking (on, on), 
neither clicking (off, off) or just one clicking (on, off) or 
(off, on). The probabilities corresponding to these events 
must sum to unity. 

The (off, off) case is identical to the case where 
both PNRDs register no photons. Defining y = 
[a 2 + 7 2 ) A, y' = (a' 2 + 7 /2 ) A, the probability in this 
case is 



P(off,off) 



1-A 2 

1 - yy' 



(42) 



The (off,on) case is also analytically tractable when in- 
corporating losses. From Eqn. (25), P (off, on) is given 

by 



P(off,on) = ^P(0,i') 



(l-A 2 )y(A-yQ 
(1 - y\) (1 - yy') 



(43) 



Note that P (on, off) is identical to the result above 
under exchange of primed and unprimed parameters. 
The probability of both detectors registering pho- 
tons P (on, on) can be tackled by using the relation 
E P ,^(on,off) = 1. yielding 



P (on, on) 



(y - A) (y' - A) (l - yy'X 2 ) 
■ (i _ yX ) (l - y'A) (1 - yy') 



Lossless subtraction 



(44) 



V. SUBTRACTION WITH THRESHOLD 
DETECTORS 

We now proceed to study the enhancement of entangle- 
ment that can be effected by subtracting photons locally 
using threshold detectors. We consider avalanche pho- 
todiodes (APDs), detectors that click when at least one 
photon is incident. Starting from Eqn. (23), the effect of 
such detectors are described by the summations over d 
from t = t max = for "off" or from t = 1 to t max = 00 
for "on", and similarly for t f . This analysis can be ex- 
tended to higher photon number threshold detectors by 
changing the ranges corresponding to the "on" and "off" 
detection. A physical example of such a detector would 



The idealised case of lossless symmetric subtrac- 
tion using APDs with unit efficiency is analytically 
tractable [16]. However, lifting the degeneracy of the 
parameters means the submatrices Ck are no longer 
centrosymmetric, rendering intractable the analytical 
method used previously. All subsequent results are there- 
fore numerical. 

In Ref. [ ] , it was shown that for squeezing parameters 
A < 0.95, the results are essentially equivalent to those 
with PNRDs. The same is true for perfect asymmetric 
subtraction: the case where (1,0) photons are simulta- 
neously detected on each mode respectively, is, within 
the limits of our numerical analysis, identical to the case 
where (on, off ) is detected with APDs, as shown in Fig. 11 
up to A = 0.8. 
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FIG. 8: Lossy symmetric subtraction with PNRDs. The (a) gain, (b) probability and (c) rate of entanglement increase as a 
function of loss 7 2 ,7 /2 , given one photon subtracted from each mode. 




FIG. 9: Lossy asymmetric subtraction with PNRDs. (a): optimised subtraction parameter cto pt and (b): squeezing parameter 
A op t under loss given asymmetric (1,0) subtraction. 



C. Lossy subtraction 

As with the PNRD case calculated above, we can use 
entanglement gain rate as the figure of merit to be op- 
timised, from which the optimal parameters can be ex- 
tracted. The optimal parameters a^ pt and A opt are shown 
in Fig. 12, with the resulting gain, probability and rate 
identical to that shown for PNRDs in Fig. 10, as a func- 
tion of the asymmetric loss parameters 7 2 , j' 2 . Note that 
the optimal parameters are also the same for both the 
PNRD and APD cases. 



1. Imperfect Asymmetric Detection 

As with perfect subtraction, the behaviour of APDs 
and PNRDs under loss is the same at experimentally ac- 
cessible values of A. In the lossless case, we have shown 
that asymmetric subtraction always produces a higher 
rate of entanglement gain, (Sec. IV C, Fig. 10). We now 
ask whether one strategy is always better than the other 
for all losses. This turns out not to be the case, as 



shown in Fig. 13. While the gain produced when sub- 
tracting from both modes is always greater than when 
subtracting from a single mode (Fig. 13(a)), for a partic- 
ular range of losses on each mode 7,7', the probability 
favours subtracting from just a single mode. This result 
is repeated in the plot showing the comparative rate of 
gain (Fig. 13(c)). 



VI. CONCLUSION 

We have investigated photon subtraction from a two- 
mode squeezed state as a means to probabilistically in- 
crease entanglement under LOCC. By defining the en- 
tanglement gain rate as our figure of merit, we are able 
to optimise over the subtraction beam splitter transmis- 
sivity T2 to maximise this quantity, in the presence of 
unequal losses on each mode. Our results may be sum- 
marised as follows: 

1. If the losses are above a threshold, which depends 
on the number of photons to be subtracted, local 
subtraction cannot enhance entanglement. 
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FIG. 10: Lossy asymmetric subtraction using PNRDs. (a) gain, (b) probability and (c) rate of entanglement gain as a function 
of the losses 7 2 , r ) 12 on modes one and two respectively. The probability of successful subtraction decreases very quickly with 
loss, and the positive gain region is bounded asymmetrically, being roughly twice as sensitive to losses on the unsubtracted 
mode. 




(a) (b) (c) 



FIG. 11: Comparing lossless subtraction strategies between APDs and PNRDs. (a) subtraction gain, (b) probability and (c) rate 
for lossless PNRDs and APDs in various subtraction strategies. The results for APDs and PNRDs for simialr strategies, up to 
A = 0.8, are identical. 



2. When it can, (1,0) subtraction seems to be the best 
strategy, whether APDs or PNRDs are used to de- 
tect the subtracted photon. 

3. However, depending on how losses are distributed 
across the modes, symmetric subtraction may be 
advantageous. 

4. Subtracting more photons produces marginal en- 
hancement in entanglement and is far less proba- 
ble. 

5. APDs are essentially equivalent to PNRDs for most 
of the A regime, including when accounting for 
asymmetric losses. 

Given any realistic scenario with lossy transmission 
channels and imperfect detectors, our conclusions out- 
line the most suitable strategy that must be adopted to 
achieve entanglement enhancement in CV systems most 
effectively. This approach also specifies the initial squeez- 
ing parameter A opt for which this particular protocol 
works best. The figure of merit we use in this paper 
is the entanglement gain rate, relevant to quantum com- 
munication applications. For other applications, it may 



be the final entanglement, or indeed the gain alone, that 
is more important. The methods presented in this paper 
can be easily modified for those purposes. 

Photon subtraction is one of the simplest operations 
introducing non-Gaussianity, thereby opening the gate 
to a large class of CV quantum information processing 
protocols such as entanglement distillation that are not 
possible in the Gaussian regime. We hope that our work 
will inform future efforts to increase entanglement via 
this technique under realistic experimental conditions. 
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FIG. 12: Lossy subtraction with APDs. (a), (b): subtraction parameters a 2 , a/ 2 , respectively, and (c): squeezing parameter A, 
optimised to produce the highest entanglement gain rate under loss given symmetric subtraction with APDs. 




FIG. 13: Comparing symmetric and asymmetric subtraction with APDs under loss. The difference between (on, on) and (on, off) 
subtraction events in (a) gain, (b) probability and (c) rate as a function of the losses 7 2 ,7 /2 on modes one and two respectively. 



[1] C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schu- 
macher, Phys. Rev. A 53, 2046 (1996). 

[2] P. G. Kwiat, S. Barraza-Lopez, A. Stefanov, and 
N. Gisin, Nature 409, 1014 (2001). 

[3] J. Eisert, S. Scheel, and M. B. Plenio, Phys. Rev. Lett. 
89, 137903 (2002). 

[4] T. Opatrny, G. Kurizki, and D.-G. Welsch, Phys. Rev. 
A 61, 032302 (2000). 

[5] L.-M. Duan, G. Giedke, J. I. Cirac, and P. Zoller, Phys. 
Rev. Lett. 84, 4002 (2000). 

[6] J. Fiurasek, L. Mista, and R. Filip, Phys. Rev. A 67, 
022304 (2003). 

[7] D. E. Browne, J. Eisert, S. Scheel, and M. B. Plenio, 

Phys. Rev. A 67, 062320 (2003). 
[8] B. Hage, A. Samblowski, J. Diguglielmo, A. Franzen, 

J. Fiurasek, and R. Schnabel, Nature Physics 4, 915 

(2008). 

[9] R. Dong, M. Lassen, J. Heersink, C. Marquardt, R. Filip, 
G. Leuchs, and U. L. Andersen, Nature Physics 4, 919 
(2008). 

[10] B. Hage, A. Samblowski, J. DiGuglielmo, J. Fiurasek, 
and R. Schnabel, Phys. Rev. Lett. 105, 230502 (2010). 

[11] R. Dong, M. Lassen, J. Heersink, C. Marquardt, R. Filip, 
G. Leuchs, and U. L. Andersen, Phys. Rev. A 82, 012312 
(2010). 

[12] H. Takahashi, J. S. Neergaard-Nielsen, M. Takeuchi, 



M. Takeoka, K. Hayasaka, A. Furusawa, and M. Sasaki, 

Nature Photonics 4, 178 (2010). 
[13] P. T. Cochrane, T. C. Ralph, and G. J. Milburn, Phys. 

Rev. A 65, 062306 (2002). 
[14] S. Olivares, M. G. A. Paris, and R. Bonifacio, Phys. Rev. 

A 67, 032314 (2003). 
[15] A. Kitagawa, M. Takeoka, M. Sasaki, and A. Chefles, 

Phys. Rev. A 73, 042310 (2006). 
[16] S. L. Zhang and P. van Loock, Phys. Rev. A 82, 062316 

(2010). 

[17] C. Navarrete-Benlloch, R. Garcia-Patron, J. H. Shapiro, 
and N. J. Cerf, Physical Review A 86, 012328 (2012). 

[18] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[19] M. Horodecki, P. Horodecki, and R. Horodecki, Physics 
Letters A 223, 1 (1996), ISSN 0375-9601. 

[20] M. B. Plenio, Phys. Rev. Lett. 95, 090503 (2005). 

[21] S. L. Braunstein and P. van Loock, Rev. Mod. Phys. 77, 
513 (2005). 

[22] M. A. Nielsen, Phys. Rev. Lett. 83, 436 (1999), URL 
http : //link . aps . org/doi/10 . 1 103/PhysRevLett . 83 . 
436. 

[23] Digital Library of Mathematical Functions. Release date: 
2012-03-23. National Institute of Standards and Technol- 
ogy from http : //dlmf . nist . gov/26 . 4. 

[24] A. Cantoni and P. Butler, Linear Algebra Appl. 13, 275 
(1976). 



14 

[25] J. Eisert, F. G. S. L. Brandao, and K. M. R. Audenaert, 
New Journal of Physics 9, 46 (2007). 



